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Abstract. We study analytically and numerically the noise-induced transition between an absorbing and 
an oscillatory state in a Duffing oscillator subject to multiplicative, Gaussian white noise. We show in a 
non-perturbative manner that a stochastic bifurcation occurs when the Lyapunov exponent of the linearised 
system becomes positive. We deduce from a simple formula for the Lyapunov exponent the phase diagram 
of the stochastic Duffing oscillator. The behaviour of physical observables, such as the oscillator's mean 
energy, is studied both close to and far from the bifurcation. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.10.Gg 
Stochastic analysis methods (Fokker- Planck, Langevin, etc.) - 05. 45. -a Nonlinear dynamics and nonlinear 
dynamical systems 
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1 Introduction 

The stability of a dynamical system can be strongly af- 
fected by the presence of uncontrolled perturbations such 
as random noise pQ . Noise-induced bifurcations are of pri- 
mary interest in many fields of applied mathematics and 
physics and occur in various experimental contexts, as for 
instance: electronic circuits 0, mechanical oscillators [3], 
laser physics surface waves 0], thermal and electrohy- 
drodynamic convection in fluids or liquid crystals .5,6,7 
and diffusion in random media 8 . In nonlinear dynami- 
cal systems, the interplay of noise and nonlinearity may 
produce unusual phenomena: noise can shift the bifurca- 
tion threshold from one phase to another 9 , it can induce 
new phase-transitions in spatially-extended systems jl(J) . 
or create spatial patterns (for a recent review see [Pi]). 

One of the simplest systems that may serve as a para- 
digm for noise-induced transitions is the nonlinear oscilla- 
tor with parametric noise |1HIU4U15| . i.e., with a frequency 
that fluctuates randomly with time around a given mean 
value [THIH?] . In a recent work [TcHTT?] . we studied the time 
asymptotic behaviour of such an oscillator in the small 
damping limit. We showed that physical observables grow 
algebraically with time until the dissipative time scale is 
reached; then the system settles into an oscillatory state 
described by a stationary probability distribution func- 
tion (PDF) of the energy E: the zero equilibrium point 
(E = 0) is unstable in the small damping limit. In the 
present work, we investigate the stability of the origin for 
arbitrary damping rate and noise strength. 



For random dynamical systems, it was recognized early 
on that various 'naive' stability criteria |16U17| . obtained 
by linearizing the dynamical equation around the origin, 
lead to ambiguous results. (This is in contrast with the 
deterministic case for which the bifurcation threshold is 
obtained without ambiguity by studying the eigenstates of 
the linearized equations [201 ■) In fact, it was conjectured 
in |17| and proved in |21| that, in a linear oscillator with 
arbitrarily small parametric noise, all moments beyond a 
certain order diverge in the long time limit. Thus, any 
criterion based on finite mean displacement, momentum, 
or energy of the linearized dynamical equation is not ade- 
quate to insure global stability: the bifurcation threshold 
of a nonlinear random dynamical system cannot be deter- 
mined simply from the moments of the linearized system. 
In practice, the transition point is usually calculated in a 
perturbative manner by using weak noise expansions jSlEl 

mm 

In the present work, we shall obtain in a non-perturba- 
tive manner the bifurcation threshold of a stochastic Duff- 
ing oscillator with multiplicative, Gaussian white noise. 
We use a technique described in ^B] where the stochas- 
tic dynamical equations are expressed in terms of energy- 
angle variables. Thanks to a detailed analysis of the asso- 
ciated Fokker-Planck equation, we shall show that the bi- 
furcation occurs precisely where the Lyapunov exponent of 
the linear oscillator changes sign. This finding is confirmed 
by numerical simulations. The stability of the fixed point 
E = of the (nonlinear) Duffing oscillator with para- 
metric noise is deduced, in this sense, from the linearized 
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dynamics. We shall also derive the scaling behaviour of 
physical observables in the vicinity of the bifurcation with 
respect to the distance to threshold. Finally, we study the 
behaviour of observables far from the bifurcation, in the 
small damping or strong noise limit. 

A rigorous mathematical theory of random dynamical 
systems has been recently developed and theorems relat- 
ing the stability of the solution of a stochastic differential 
equation to the sign of Lyapunov exponents have been 
proved |22|. In the present work, we do not rely upon 
the sophisticated tools of the general mathematical the- 
ory but follow an intuitive approach, easily accessible to 
physicists, based upon a factorisation hypothesis of the 
stationary PDF in the limit of small and large energies. 

This article is organized as follows. In Section we 
define the system considered, present the relevant phe- 
nomenology, and introduce energy-angle variables for the 
stochastic Duffing oscillator. We study the stationary PDF 
of the energy in Section [3J In Section 0] we derive the sta- 
bility criterion for the stochastic Duffing oscillator and 
obtain its phase diagram. In Sectional we perform a local 
analysis of the physical observables in the vicinity of and 
far from the bifurcation threshold. Sectionals devoted to 
some concluding remarks. In Appendix^ we prove some 
useful identities satisfied by the Lyapunov exponent. In 
Appendix El an explicit formula for the Lyapunov expo- 
nent is derived. 



2 The Duffing oscillator with multiplicative 
white noise 

2.1 Notations and phenomenology 

A nonlinear oscillator with a randomly varying frequency 
due to external noise can be described by the following 
equation 

d 2 x dx , o _ .. dU(x) 

where x(t) is the position of the oscillator at time t, U(x) 
an anharmonic confining potential and 7 the (positive) 
friction coefficient. The linear frequency of the oscillator 
has a mean value u> and its fluctuations are modeled by 
a Gaussian white noise of zero mean- value and of 
amplitude T> 

(£(*)> = 0, 

(atw)) = vs(t-t>). (2) 

In this work, all stochastic differential equations are inter- 
preted according to Stratonovich calculus. We shall study 
the Duffing oscillator with multiplicative noise, the con- 
fining potential IA being given by 



(a) 
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(b) 
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(c) 
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Fig. 1. Phase plane plots of a typical trajectory of the noisy 
Duffing oscillator. Eq. is integrated numerically with a time 
step 5t = 5 I0" 4 . (a) a = 0.5, A = 20, 10 < t < 100; (b) a = 
2.0, A = 20, t < 100; (c) a = 2.0, A = 60, 10 < t < 100. Early- 
time relaxation towards noisy oscillations has been omitted for 
clarity in cases (a) and (c). The trajectory spirals down towards 
the globally stable origin in case (b). 



U(x) = A — 
4 



(3) 
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By rewriting time and amplitude in dimensionless units, 
t := tut and x := \ 1 ^ 2 L)~ 1 x, respectively, Eq. becomes 



d 2 x dx a 



(4) 



where El (t) is a delta-correlated Gaussian variable: 

<-(*)> =0, 
(S(t)S(t'))=AS(t-t'). (5) 



The parameters 

a = y/u> and A = T)/uj 3 



(0) 



correspond to dimensionless dissipation rate and to noise 
strength, respectively. In the rest of this work, all physical 
quantities are expressed in dimensionless units. 

In Fig. ^ we present, for three values of the parame- 
ters, a trajectory in the phase plane (x, v = x) character- 
istic of the oscillator's behaviour. Here as well as for other 
simulations presented in this article, Eq. |0J is integrated 
numerically by using the one step collocation method ad- 
vocated in and described in detail in [IB]. Initial con- 
ditions are chosen far from the origin, with an amplitude 
of order 1. Noisy oscillations are observed for small values 
of the damping parameter a [See Fig. Q]a]. A larger value 
of a turns the origin into a global attractor for the dy- 
namics [See Fig.[J]b]. Increasing the noise amplitude A at 
constant a makes the origin unstable again [See Fig.r^c]. 

For all values of the parameters (a, A), the origin (x — 
x = 0) is a fixed point of the random dynamical sys- 
tem (0J. Our main goal, in this work, is to investigate the 
stability properties of the fixed point at the origin when 
the control parameters vary. 



2.2 Energy-angle representation of the Duffing 
oscillator 

We now follow the method presented in ^H] and introduce 
the energy-angle variables associated with the Duffing os- 
cillator. The deterministic and undamped oscillator cor- 
responding to Eq. is an integrable dynamical system 
because the mechanical energy E, defined as 



E=-x 2 



1 2 
— X 

2 



1 4 

r 



(7) 



is a conserved quantity. Introducing the angular variable 
4>, the constant energy orbits of the Hamiltonian system 



are parametrized as 



4E 2 



1/4 



AE+ 1 

1/2 cn(</> ;fc) 



sd (4> ; fc) 



(225) 



dn (</> ; fc) ' 



(8) 
(9) 



where the functions sd, cn and dn are Jacobi elliptic func- 
tions [22] • The elliptic modulus fc is a function of the en- 
ergy and is given by 



V4£ +1-1 



2V4£ + 1 



(10) 



In the (E, <j)) coordinate system, the dynamical equations 
for the Duffing oscillator, without noise and dissipation, 
become simply 



E = 0, 

<P = (4E + 1)< 



(11) 



We reintroduce the multiplicative noise and the dissi- 
pation terms, and obtain a system of coupled stochastic 
equations that describes the evolution of the energy and 
angle variables: 

E = J E {E,4>)+D E {E,4>)E{t), (12) 
<P = J^E,4>)+D^(E,4>) S(t) , (13) 

where the drift and diffusion coefficients J E , J^, De and 
£>0 are defined by (see rTS| for details of the derivation) : 



J E (E,ct>) = -2aE^Ml 



D E (25, i 



dn (tfi; fc) 
2E sn(</>; fc)cn(</>; fc) 



(4£+l)3 dn 3 (0;fc) 



(14) 
(15) 



. ,. , ,i 2E + 1 sn(d>: k)cn(d>: fc) . 

J,( B ,»)-(« B + l ) , +a — g^L-gi 

E 4 aE cm>;fc)sd 3 (<£;fc) 

-sd (0; fc) + 



(4E+1 



(4E + l)t dn 2 (</>;fc) 



+ 



E sn(4>; fc)cn(</>; fc) 



(4E + 1)I dn 3 (<?(>; fc) 7 
q£ cn 2 (0;fc) 



sd 2 (6<; k)dS 



(4£+l)t dn 4 (0;fc) Jo 



s<T(0;A;)ef0. 



(17) 



3 Analysis of the Fokker-Planck equation 

Using the same notations, the Fokker-Planck equation that 
governs the time evolution of the probability distribution 
function Pt(E,(f>) associated with the stochastic differen- 
tial system ([T21)-([T3l reads 25 



d t Pt(E,< 



-Be (JeP) - 84, (J+P) 



(18) 



+ y [d E (Z) E d E (D B P))+d E (D E d (D P)) 
+9 {D^d E (D E P)) + (D+diiD+P)) } . 

The goal of this section is to determine the conditions 
under which there exist stationary solutions P s tat(2?, (f>) of 
Eq. 1|18J) with a non-zero mean energy. Such an extended 
PDF describes an oscillatory asymptotic state. 

Obtaining an analytical expression of P sta t(P, </>) in the 
general case is a daunting task. In the following, we shall 
perform a local analysis of Eq. I|18fl for small and large 
values of E. In these limiting cases, the expressions of the 
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drift and diffusion coefficients take simpler forms. We may 
write, in full generality, 



Pstat 



Pstat (</>|P) Pstat (P), 



(19) 



where P s tat(<^|P) represents, in the stationary state, the 
conditional probability of the angle ^ at a fixed value 
of the energy P. Assuming the conditional probability 
Pstat {4>\E) to be independent of P when E<1 and E ^> 1 
(respectively Sees. I3.ll and 13. 2|) . we shall derive the local 
behaviour of a stationary solution of the Fokker-Planck 
equation i|18|) . This solution is a legitimate PDF when it 
is normalizable: the normalizability condition provides the 
location of the bifurcation 1261. 



3.1 Small E limit: the linearized stochastic oscillator 

When the oscillator's mechanical energy E is small, the 
position x(t) and the velocity x{t) are also small: nonlinear 
terms may be neglected. The Duffing oscillator simplifies 
to a linear (harmonic) oscillator: 



x + ax + x = x 5 (t) 



(20) 



with (~(t)S(t')) = A S(t - t'). The energy now reads 
E = \ x 2 + | x 2 , and Eqs. (JSJ) and I© reduce to 



x = v 2E sin (, 
x = V2P cos ( 



(21) 
(22) 



In the small E limit, the elliptic modulus k goes to zero, 
the elliptic functions reduce to circular functions [2J , and 
the drift and diffusion coefficients simplify to 



Je{E,4>) = -2aE cos 2 (f> 
D E (E,(j)) — 2E sin</>cos0, 
J<f,(E,(f)) = 1 + a sin 4> cos 4> > 
D<i,(E,<t>) = -sm 2 (f>. 



(23) 
(24) 
(25) 
(26) 



Substituting the expressions lj2l5|) - (|2l)|) in Eq. (|T%|) and in- 
tegrating over P, we obtain an independent Fokker-Planck 
equation for the marginal PDF Pt{4>): 



d t Pt{(f) = -d<t, (1 + a sin (/> cos 0) P t {(j>) 



A 



sm 2 <^ (sin 2 <^>P t (</,)) . (27) 



An explicit formula for the stationary angular measure 
-Pstat (</>), solution of Eq. (|27|l . is derived in Appendix iBl 

We now make the hypothesis that in the small E limit, 
the stationary conditional measure P s t a t(0|P) becomes in- 
dependent of E and is identical to the stationary angular 
distribution P s tat(0) of the harmonic oscillator: 

for P«l, P stat (0|P) =P sta t(0). (28) 
Equation (|19|1 becomes 

Pstat {E,</>) = Pstat (P)PstatW>), (29) 



i.e., we assume that the stationary PDF factorizes when 
E <C 1. An exactly solvable equation for P s t a t(E) can now 
be obtained by inserting expressions (|23fl - l|26[l and the 
identity 1)29(1 in the Fokker-Planck equation l(18|l and then 



averaging over 



Y (sin 2 <j> cos 2 • 



d E (EP stat {E) 



(30) 



(4 (sin 2 (sin 2 <\> - cos 2 <£)) - a (cos 2 <f>)) P stat (P) . 

Expectation values denoted by brackets ( ) are calculated 
using the stationary angular measure P s tat(0)- Introduc- 
ing the notations: 

fj,(a,A) = y (sin 2 cos 2 0) , (31) 

A(a, A) = y (sin 2 (sin 2 4> - cos 2 (f>)) - a (cos 2 </>)(32) 

Eq. <|30|l becomes: 

[i d E (pP sta t(P)) =A P stat (P), (33) 

and admits the solution: 

P stat (P) cx P^ 1 for P«l. (34) 

Thus, for E <C 1, the PDF P s tat(P) behaves as a power 
law of the energy. 

The coefficient fi is always positive, but the sign of A is 
a function of the parameters a and A. For negative values 
of A, the stationary solution 1341) is not normalizable and 
cannot represent a PDF; the only stationary PDF solution 
of (|18J) is then a delta function centered at the origin: 



Ptat(P) =S(E) for yl < 0. 



(35) 



Thus for negative values of A, the origin is a global at- 
tractor for the stochastic nonlinear oscillator. Non-trivial 
stationary solutions exist only for A > 0: the bifurca- 
tion threshold corresponds to values of the parameters for 
which A(a, A) = 0. 



3.2 Large P limit 

For very large values of P, the elliptic modulus fc, denned 
in Eq. (|10(l . tends to l/\/2, and the drift and diffusion 
coefficients become 



Je(E,< 



-2aE 



cn 2 (</>) 
dn 4 (0) ' 



(36) 
(37) 



p £ (p,0)^V2P^^M, 

dn (4>) 

J (P, 0) = (4P) ' + - dn 3^^ ■ ( 38 ) 

J^(£,0) = - ^ i/4 sd 2 (0), (39) 
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in agreement with Eqs. (3) and (4) of JHj (we have omit- 
ted the elliptic modulus 1 /y/2 for sake of simplicity) . Using 
expressions <(36() to (|39f) for the drift and diffusion coeffi- 
cients, the stochastic equations ((12(1 and (|13|) reduce, in 
the large E limit, to 



E = -2aE 



{AEf 



dn 4 (0) 
/4 ^ a sn(0) cn{4>) 
2 dn 3 (0) 



dn J (</>) 



In terms of the variable Q — E 1 / 4 , Eq. 140(1 becomes linear 
and has a form similar to Ornstein-Uhlenbeck's equation. 
Because the elliptic functions are bounded, fl saturates to 
a value of the order of ~. Thus, we deduce from Eq. ijlTJl 
that the phase <f> grows linearly with time. Hence, is a 
fast variable as compared to J?. Integrating out this fast 
variable from the Fokker-Planck equation leads to an ef- 
fective stochastic dynamics for i?, as explained in detail 
in |I8j . An equivalent formulation is to assume that, in 
the large E limit, the stationary conditional probability 
-Pstatt^l^), defined in Eq. (|f 9fl . becomes uniform in <j> and 
independent of E, i.e., 



for E > 1 , P si 



\E)~1 



(42) 



Averaging the dynamics over <f> yields the following ex- 
pression for the stationary PDF |13U19| : 



for E » 1, P stat (£) cx E-^cxp j 



8 a E 1 / 2 
3MA 



(43) 



with M = (3/4) 4 ~ 0.3655, T being the Eulerian 

factorial function. A stationary solution of l|18(l decays as 
a stretched exponential and is always integrable at infinity 
irrespective of the values of a and A. 



4 Phase Diagram of the Duffing oscillator 

In the long time limit, the energy of a harmonic oscilla- 
tor with parametric noise (121 )[) varies exponentially with 
time 16 . This behaviour is characterized by a Lyapunov 
exponent, defined as lim^oo ^Togi?, where the over- 
line means averaging over Pt(E, <fi) 27 . [The factor 1/2 is 
included for consistency with the usual definition of the 
Lyapunov exponent in terms of the position x(t).] 

In Appendix lAl we prove that the coefficient A(a, A) 
defined in Eq. 1(32(1 is equal to the Lyapunov exponent: 



f 



A(a, A) = lim — log E . 

t— too It 



(44) 



For the (nonlinear) Duffing oscillator, we found in Sec- 
tion |3 that a non-trivial stationary solution P s t a t(E, </>) of 
Eq. I(I8() exists only when the Lyapunov exponent A(a, A) 
of the associated linear oscillator is positive. This extended 
PDF describes an oscillatory asymptotic state for the Duff- 
ing oscillator, as opposed to the absorbing state associated 



100 r 
80 
60 



A > 

Oscillatory state 




40- 



20- 



Fig. 2. Phase diagram of the Duffing oscillator with multi- 
plicative noise. The critical curve separates an absorbing and 
an oscillatory phase in the (a, A) plane. The dashed line rep- 
resents the asymptotic approximation, Eq. 147H . 



with the delta function 1(35(1 . The location of the stochastic 
bifurcation threshold of the Duffing oscillator, which sep- 
arates the two states, is given by the curve A(a, A) = 0. 
The following explicit formula is derived in Appendix iBl 



A(a,A) = 



1 J Jo °° du 



>u e 



r+°° Au „-- 
Jo ,AT 



(1- 



-)"+ Ta 



(45) 

The transition line is determined by the relation be- 
tween a and A for which the Lyapunov exponent vanishes. 
For a given value of A, the critical value of the damping 
a = a c {A) satisfies 



/ o +oc d M 0Iexp{-i((l-4)u+f;)} 
/o +0 °^xp{-i((I-^) U+ f!)} ' 



(46) 



The critical curve a — a c (A) [Eq. 1(46(1 ] is represented in 
Fig. [3 It separates two regions in parameter space: for 
a < a c (resp. a > a c ) the Lyapunov exponent is positive 
(resp. negative), the stationary PDF is an extended func- 
tion (resp. a delta distribution) of the energy, the origin 
is unstable (resp. stable). For a detailed analytical study, 
we must now distinguish three cases corresponding to an 
undcrdamped (a c < 2), critically damped (a c — 2) and 
overdamped (a c > 2) oscillator. 



4.1 Underdamped oscillator 



The case a c < 2 can formally be mapped onto that of a 
linear oscillator with no damping. In terms of the variables 



1- ^ 



-1/2 



and Z\<- = A 



1 -3/2 



Eq. gSJ 
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1- — 



A< 



i/ 2 J + -d,>exp{-^(, + ^)} 
Er=o^^(3-+|) (#)' 



oo (-1)* 



n=0 



r(3n 



2/ ^48 J 



Retaining only the lowest order term in this series, we 
obtain the limiting behaviour of the critical curve a c (A), 
for small values of a and A: 



A 



(47) 



The change of variables used here breaks down at a c — 2, 
we shall therefore consider this case separately. 



4.2 Critically damped oscillator 

For a critically damped oscillator, we substitute a c — 2 in 
Eq. and obtain 

| + °°^vWp(-^) r (i) 

2 r+-du_f_ J!L \ ~ {6Ac) TU)- (48) 

Jo 7s ex P^ 6Z7j Ve; 



Thus, the critical noise amplitude A c for a c — 2 is 
^ = 3^375- - 41.2969..., 



(49) 



and the expression i|45|) for the Lyapunov exponent re- 
duces to 

Z\\ 1/3 

(50) 



A = { A~ ] ' 



4.3 Overdamped oscillator 

For a c > 2, we define v = 
r 2 n-3/2 

A % - 1 



-1/2 



and Z\- 



0£ 

4 



Equation (|46|) then becomes 
1/2 _ f Q +c °dv^e-^(- v+ £) 



r +oo 



-D+4 



n=0 n! 1 V 6 7 1^3^ ^ 



(51) 

n/3 



(64>) 



Y^oo j_ r (2n±l\ ( 48 \ 
Z^n=0 n! V 6 / \ J 



n/3 



When a c — > oo, the left hand side of Eq. I|51|) the pre- 
vious equation converges to 2 and Z\> tends to a limiting 
value K, that we computed numerically: K, ~ 28.337. Thus, 



in the large a c limit, we find the following relation between 
a c and A 



A~K 



-i 3/2 



^a? ~ 3.54 a? 



(52) 



When a c — 2 + , the left hand side of Eq. I|51|) increases 
to infinity. Retaining only the leading order term of the 
series, we obtain 



-1/2 



This equation is consistent with Eq. (|49|l . 



5 Physical observables 

The goal of this section is to determine how physical ob- 
servables, such as the energy and the position of the os- 
cillator, scale near the bifurcation threshold as well as far 
from it. The agreement between analytical calculations 
and numerical simulations will validate the assumptions 
made in Section [3] concerning the behaviour of the condi- 
tional PDF P s t a t((f>\E) for small and large values of E. 

Let us first of all review the results obtained in Sec- 
tion |3| For parameter values such that A < 0, we demon- 
strated that Pstat(E) = 5(E): the origin is a global attrac- 
tor of the dynamics. Conversely, when A > 0, the origin 
becomes unstable. The stationary measure of the energy 
is then an extended function: it behaves as a power law 
for small values of E [see Eq. I|34|) ] . and decreases as a 
stretched exponential for large values of E [see Eq. I|43|) ]: 



for P< at (E) = AE^ 1 , 

for £»1, P> at (E) = BE-iexp 



" 3MA 



(54) 



The precise values of A and B can only be obtained by 
solving the full Fokker-Planck equation (|18fl . and match- 
ing the two asymptotic expressions in the intermediate 
regime E ~ 1. 



5.1 Behaviour in the vicinity of the bifurcation 

For a given value of the damping rate a, let us choose the 
noise strength A such that < A — A c (a) <C 1. In this 
case, the parameters of the system are tuned just above 
the bifurcation threshold and the Lyapunov exponent is 
slightly positive: < A -C 1. The stationary PDF is an 
extended function but, as can be inferred from Eq. 
most of its mass is concentrated in the vicinity of E = 0. 
Therefore, a good approximation of the stationary PDF 
is obtained by taking it equal to P^ at (E) over a finite 
interval in the vicinity of (e.g., < E < Eq), and equal 
to zero outside this interval. The mean value of the energy 
is given by: 



(Eh 



rEo pi 
Jo E " 



dE A 

oc E Q — 

' x dE M 



(55) 
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Fig. 3. Scaling behaviour of physical observables of the 
stochastic Duffing oscillator in the vicinity of the bifurca- 
tion threshold as a function of the reduced control parameter 
e — (A — A c )/A c . (a) Average energy; (b) Variance of the po- 
sition. Dashed lines corresponding to a linear function of e are 
drawn to guide the eye. The corresponding linear-linear graphs 
are shown as insets. Eq. @ is integrated numerically for a — 2, 
with a time step St — 5 10~ 4 . The threshold value is given by 
Eq. l|3nj: A c (2) ~ 41.3. Ensemble averages are computed over 
10 3 realisations. 



If we now take into account that: (i) A(a, A c (a)) = 0; 
and (ii) the Lyapunov exponent and the coefficient \i are 
regular functions of the noise strength A when A ^ 0, we 
obtain from Eq. (|55|) 
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Fig. 4. Scaling behaviour of the mean energy of the general 
nonlinear oscillator in the vicinity of the bifurcation thresh- 
old as a function of the reduced control parameter e = (A — 
A c ) /A c . Dashed lines indicating a linear dependence in e are 
drawn to guide the eye. Eq. Q is integrated numerically for 
a = 2, and with U(x) = x 2 "/(2n), n = 2, 3 and 4. Data ob- 
tained for n = 2 are the same as in Fig. |3] The time step is 
St — 5 10 -4 ; ensemble averages are computed over 10 3 realisa- 
tions. 



we find 



<* 2 > S 



J f -Pstat (E, 0) 2E sin 2 (j, dE d<P 



= 2 (sin 2 4>) (E) sUt 
on (A- A c ) . 



(57) 



The variance of the position scales linearly with the dis- 
tance to threshold [See Fig.|3|b] . Using a similar argument, 
one also finds that (u 2 ) s tat oc (A — A c ). 

As seen in Fig. 01 the linear scaling (|56|l of the energy 
with respect to the distance to threshold is observed for 
U(x) oc a; 2 ™, n — 2, 3 and 4, where A c (a) is always given 
by A(a, A c (a)) = 0. The threshold is independent of the 
nonlinearity involved in the confining potential hi appear- 
ing in Eq. 0), because the local expression P^ at (E) of 
the stationary PDF in the vicinity of 73 = is obtained 
from the linearized stochastic differential equation. The 
observed dependence of the constant of proportionality in 
Eq. H56f) upon the stiffness of the anharmonic potential 
remains to be understood. 



(E) stat oc A{a, A) - A{a, A c ) oc{A- A c ) , (56) 



where the last identity is obtained by a Taylor expansion 
of A for A close to A c . This result is confirmed by numer- 
ical simulations [See Fig. El a]. A similar reasoning shows 
that (£%tat scales linearly with a — a c (A) at fixed noise 
strength A. 

Let us now consider the variance of the oscillator's po- 
sition close to the bifurcation. Using Eqs. Ij21|l and (|29|) . 



5.2 Scaling far from the bifurcation 

We shall now study the scaling of the physical observables 
far from the bifurcation point. This situation is encoun- 
tered, for instance, when the damping rate is vanishingly 
small for a given noise strength. The Lyapunov exponent 
is then a finite, strictly positive, number. When a — ► 0, 
most of the mass of the stationary PDF is concentrated 
at large values of E as can be seen from Eq. ()54(l . Thus, 
a satisfactory approximation of the stationary PDF is to 
assume it to be identical to the expression P s ^ at (E) for all 
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Fig. 5. Scaling behaviour of the mean energy of the gen- 
eral nonlinear oscillator far from the bifurcation as a func- 
tion of the parameter A /a. Eq. Q is integrated numerically 
with U(x) = x 2n /(2n), n = 2, 3 and 4. Dashed lines corre- 
spond to the predicted scaling l|6()|l . with a scaling exponent 
n/(n - 1) = 2, 3/2 and 4/3. The time step is St = 5 10" 4 ; 
ensemble averages are computed over 10 3 realisations. 



E. After a suitable normalisation, we find 



1 



2Er(3/2) V 3MA 



exp 



SaE 1 ' 2 
' 3M A 



(58) 

with the numerical factor M. defined in Eq. I|43l) . In this 
low damping regime, we recover, as expected, the results 
derived in Refs. I18lll9j . The average energy is readily ob- 
tained from Eq. (I5S 



(Eh 



135 
256' 



M 2 



0.070 



A 



(59) 



This result is at variance with the usual equilibrium fluc- 
tuation-dissipation theorem ((E) cq cx A/ a); yet this is 
not a contradiction since Eq. Q models a system out of 
equilibrium. 

Far from the bifurcation threshold, the scaling behaviour 
depends strongly on the stiffness of the potential U. If we 
take U(x) oc x 2n , Eq. 159fl can be generalized to fT!|| 



(E) s 



1 

2n 



r ( _3n+l_\ 
1 \2(n-l)J 



r 



2(n-l) 



(60) 



(n - 



r(-3-\ r ( 3n+l 
^ \ 1 \2n) 1 V 2n 

r ' ' ' 



1 \ p I 3k+3 > 
In) V 2n ' 



mechanical energy is totally dissipated in the long time 
limit and the system converges to the origin. This pic- 
ture is no longer valid once the oscillator is subject to 
parametric noise. Indeed, the dissipative Duffing oscilla- 
tor with parametric white noise undergoes a noise- induced 
bifurcation between: (?) an absorbing state reminiscent of 
the deterministic system; and (ii) an oscillatory state, re- 
sulting from a dynamical balance between energy dissipa- 
tion and (noise-induced) energy injection. Noisy oscilla- 
tions appear when the fixed point at the origin becomes 
unstable: this transition is thus referred to as a (forward) 
stochastic Hopf bifurcation. Both states may be charac- 
terized by the stationary probability distribution function 
of the oscillator's mechanical energy E. The PDF of the 
absorbing state is a delta-function, whereas the oscillatory 
state corresponds to an extended PDF. 

Our main result concerns the location of the transition 
point: a noise-induced bifurcation for the Duffing oscilla- 
tor occurs precisely where the Lyapunov exponent of the 
linear oscillator with parametric noise changes sign. We 
emphasize that this result is genuinely non-perturbative. 
It results from a simple hypothesis on the behaviour of 
the conditional probability of the angular variable at fixed 
energy, valid for very small values of the energy and for 
arbitrary noise amplitude. This assumption is physically 
plausible (i.e. mean-field-like) and justified by the conse- 
quences deduced from it. Having identified the Lyapunov 
exponent of the stochastic linear oscillator as the impor- 
tant physical quantity, we give a simple analytical expres- 
sion of this exponent as a function of the system's pa- 
rameters, the damping rate and the noise strength. The 
location of the transition line that separates the absorbing 
and oscillatory states of the (nonlinear) Duffing oscillator 
follows immediately. 

Although we have considered here only a cubic oscilla- 
tor with parametric Gaussian white noise, we believe that 
the relation between the sign of the Lyapunov exponent 
and the location of the noise-induced transition holds irre- 
spective of the nonlinearity and of the nature of the noise. 
Finally, we demonstrate that all physical observables scale 
linearly with respect to the distance to threshold. Indeed, 
the argument we use to evaluate the mean energy (E) may 
be easily generalized to show that all moments (E k ) also 
scale linearly, suggesting that the concept of critical ex- 
ponents is not relevant here. One would certainly like to 
test experimentally whether noise- induced transitions may 
indeed be characterized by a multifractal order parame- 
ter, at variance with the theory of equilibrium continuous 
transitions. 



This prediction is again confirmed by numerical simula- 
tions [See Fig. 



A Some useful identities for the Lyapunov 
exponent 



6 Conclusion 

For a deterministic Duffing oscillator with linear damping, 
the fixed point at the origin (x = x = 0) is a global at- 
tractor: irrespective of the choice of initial conditions, the 



In this appendix, we prove that the Lyapunov exponent, 
defined as lim^oo ^ log-E, satisfies the relation iflljl . Us- 
ing the 'radial' variable p defined by 



1 



log(2£) , 



(61) 
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we rewrite Eq. ((20(1 as a system of nonlinear coupled stochas- 
tic equations 

p = —a cos 2 <j) + sin <f> cos 4> S(t) , (62) 
4> = 1 + a sin (j) cos <j> — sin 2 4>E{t) . (63) 

The Fokker-Planck evolution equation for the joint prob- 
ability distribution P t (p, (f>) is then given by [2E| 

dtPt = ol cos 2 4> dpPt — 9^ [(1 + a sin <f> cos 4>) P t ] (64) 
+ ^{ sin 2 4> cos 2 (j> d 2 pp P t - sin <j> cos </> d 2 pif) [sin 2 P t ] 
[sin 3 0cos0d p P t ] +9 [sin 2 09 (sin 2 0P t )] } 

From Eq. i|64|) we derive the time evolution of the mean 
value 7j, where the overline indicates averaging with re- 
spect to Pt(p,(j)). After suitable integration by parts, we 
obtain 

^-p = —a cos 2 + ^sin 2 (sin 2 — cos 2 <p) . (65) 



B Calculation of the Lyapunov exponent 

Expressions for the Lyapunov exponent have been avail- 
able for a long time in the literature but they are 
rather unwieldy. In a recent article about the quantum lo- 
calisation problem |28| . an expression for the localisation 
length is derived: it is shown to be identical to the inverse 
Lyapunov exponent of the linear oscillator with paramet- 
ric noise and zero damping. In this appendix, we generalize 
the calculation of [2Ej to the dissipative case and derive a 
particularly simple formula for the Lyapunov exponent. 

Since A = (cotan</>) [See Eq. 1)71(1 ]. we introduce the 
auxiliary variable z defined as 



z = x/x = cotan . 
Then, equation ((20(1 becomes 

z + l+ az + z 2 = S(t) 



(72) 



(73) 



Similarly, we deduce the time evolution of log | sin < 



—log | sin <p\ = cotan <fi + a cos 2 . 
A 



— — sin 2 (sin 2 — cos 2 0) . (66) 



From Eqs. and iffi ty . we obtain 



-p = cotan i 



dog | sin< 



(67) 



dt dt 
At large times, the PDF Pt(E,<fi) reaches a stationary 



The stationary Fokker-Planck equation associated with 
Eq. J73 is given by 

-r ^-P a tat = -(l + az + z 2 )P stat + J, (74) 
2 dz 

where we have introduced the stationary current J. The 
solution of this equation is found by the method of the 
variation of constants: 

P«*(z) = ^ f exp j|(<%) - *(*))} dy, (75) 

where <P{y) = y+%y 2 + ±y 3 . Writing P stat (z) dz = 1, 
we obtain the normalisation constant N: 



limit P s t&t(E, (/)). Thus, when t — + oo, the averages log | sin0 
and cotan have finite limits, given by (log|sin0|) and 
(cotan 0), respectively, where the expectation values de- 
noted by brackets () are calculated by using the marginal 
PDF, P sta t(# Therefore 



N = 



f + °° f z (2 

/ dz dyexpl - #( z )) 

J — oc J — oo L ^ 



(76) 



Following 28 , we change the y variable to u = z — y in 
Eq. (75)1 : 



lim —p— (cotan < 
t^oo dt x 

From Eqs. (BJJ and JBjJ, we derive 



(68) 



N 



dz / dwe - ^ 
■oo Jo 



1 



1 d. 



(77) 

The integral over z is now Gaussian, and can be evaluated 
before the integral over u. Hence, we obtain 



lim — \ogE — lim — — logP = lim —p = (cotan 0) . 
t^oo 2t t^oo 2 dt t^oo dt 

(69) 

Similarly, we take the stationary limit of Eq. I(66|) and 
obtain 

(cotan 0) = —a (cos 2 <f>) + — (sin 2 (f> (sin 2 4> — cos 2 4>)) . 

(70) 

The right hand side of this equation is nothing but A [See 
Eq. (32H- This concludes the derivation of Eq. (gU : 



N 



ttA 

2 J v a 



du f 2 / a" ii 
exp <j — - [ (1 — )u + 



A V 4 

The use of the identity 

-Pstat(0)d0 = P sta t(z)dz 

yields the stationary angular measure 
1 



(78) 
(79) 



lim — logP = (cotan 0) = A(ot,A) . 



(71) 



p s tM = 



N sin 



dy e~ 



(#(y)—<f (cotan (ji)) 



(80) 
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We deduce that the Lyapunov exponent satifies 

A= I P sta t(0)cotan0d0 
Jo 



(81) 



1 

N 



+ OC 



dzz\ exp<{ — (*(y) -*(«)) }>dy. 



We again change the variables from (y, z) to (u, z) with 
u = z — y, and obtain 



1 / >+0 ° 
A = — dzz 
N 



p/ rv 2 u ^ / 2 \ 2 "\ 

due ^ 3 / 

(82) 

Evaluating the Gaussian integral in z leads to the expres- 
sion of the Lyapunov exponent given in the text 

Jo 

.'(83) 

This formula also appears in recent mathematical litera- 
ture It constitutes an exact result, obtained without 
approximations. We have computed numerically the Lya- 
punov exponent by measuring logi? in the time-asympto- 
tic regime and thus confirmed the validity of i|83|) . 
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